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A closed system of the equations for the local Bloch vectors and spin correlation functions is 
obtained by decomplexification of the Liouville-von Neumann equation for 4 magnetic particles 
with the exchange interaction, that takes place in an arbitrary time-dependent external magnetic 
' field. The numerical analysis of the quantum thermodynamic variables is carried out depending 

on separable mixed initial state and the magnetic field modulation. Under unitary evolution local 
j ■ dynamics of power production in the finite environment is investigated. 
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I. INTRODUCTION 



The classical thermodynamic heat engine converts heat energy in a mechanical work with the help classical me- 
chanical system in which gas extends and pushes the piston in the cylinder. Such heat engine receives energy from the 
high-temperature reservoir. The energy part from this reservoir is converted in mechanical work, and the part of it 
\ is transferred in the low-temperature reservoir. The classical heat engine reaches its peak efficiency if it is reversible. 

a Because of the impossibility to construct the heat engine which is ideally reversible, in 1824 Carnot [l[ has offered 
I \ the mathematical model of the ideal heat engine which is not only reversible, but is also cyclic. In the last decades 
'"^ . great efforts has been devoted to the investigation of the quantum properties of the working substance, the search 
^ ' and practical implementation of the quantum analogue of the Carnot cycle in microsystems. 

O . The operation of quantum heat engines that employ as working agents multi-level systems, for example harmonic 
oscillators, free particles in a box, three-level atoms or electrons submitted to magnetic fields, were introduced in [2[, 
Q, The two- level quantum systems similar to particles with spin 1/2 are the essential ingredients for quantum 
CN ' computation, but the coupled spin systems can also be used as quantum thermodynamic engines [l-d, [l^, El . The 
^ \ quantum analogue of the Carnot cycle requires a dynamical description of the working medium, the power output and 
the heat transport mechanism. The spin system " working gas " in an external field has its own physical properties. 
These properties at weak interaction with the environment (heat baths) can, as a rule, slightly be deformed. The 
, purpose of this paper is to predict these properties to know what we can expect at the contact with the environment 
which is usually considered in the Markovian approximation for a derivation of the equation for the reduced matrix. 

The paper is organized as follows. In Sec. II we introduce the model Hamiltonian. In Sec. Ill the Liouville-von 
Neumann equation for a density matrix for four particles with spin 1/2 with exchange interaction in variable magnetic 
' field wc write down in the Bloch representation in terms of the local Bloch vectors and spin correlation functions. We 
• • ■ describe the conservation laws which effectively supervise the numerical calculations. In Sec. IV wc describe the local 
quantum thermodynamic parameters of the spin subsystems. Our numerical results are detailed in Sec. V. For the 
separable mixed initial state we numerically investigate the quantum thermodynamics of a particle in the environment 
$^ ■ of three others depending on modulation of an operating field and initial conditions. Our results are summarized and 
5^ , discussed in Sec. VI. It will be numerically found, that under unitary dynamics the work production of one part of 
system is compensated by absorption of work produced by the other part. The work production by the subsystem is 
accompanied the entropy growth and vice-versa, the entropy becomes less with work absorption by the subsystem. 
Some necessary auxiliary details for numerical results are represented in the Appendix. 
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II. MODEL HAMILTONIAN 



The Hamiltonian of four coupled particles e,p, n, u with spin 1/2 in the external ac magnetic field h = (hi, /12, ft-a) 
looks like 

+J™<s^ + Jf"<s^ + J""sf s^), (1) 

where h^,h^,h2,hf are the Cartesian components of the external magnetic field in the energy units, operating 
on corresponding particle; = ^a^ 1^ ao (gi ctq ctq, = ^ctq (8) di (g) CTq ® CTq, s" = ^ctq (8) CTq (81 CTi CTo, 
s" = ^(To (8) (To <8) fo (8 (Ti are the matrix representation of spin operators ([T|); the Pauli matrices are equal to 

do = (^01)' ~ (10)' ^ 0*)' ''^^ ~ (o-l)' symbol of direct product [H] ; 

jep jen jpn ^ jeu ^ jpu ^ jnu g^j,g ^j^g constants of isotropic cxchangc interaction between spins; the summation over 
e,p, n, u is absent. 



III. DECOMPLEXIFICATION OF THE LIOUVILLE-VON NEUMANN EQUATION 

The Liouville-von Neumann equation for the density matrix p, describing the dynamics of four-spin system, looks 
like 

idtp^[H,p], p{t^Q) = p^. (2) 
Let us present the solution of the equation ^ as 

1 , 
P = -^Rap-iS^a (81 cr/3 (8) cr^ (8) (Jg, p ^ p, Tr p = 1, i?oooo = 1, (3) 

Hereinafter summation is taken over the repeating Greek indices from zero up to four, and over Latin indices from 
one up to three. The four coherence vectors (the Bloch vectors) widely used in the magnetic resonance theory, are 
written as 

Riom = Tr pai (8 CTo (8 (To (8 (To , (4a) 



^oioo = Tr puQ ®ai®(To®ao, 



^QOiO = Tr pcTQ (8 (To (8 (Ti (8 (To , 



(4b) 
(4c) 



Romi = Tr p(To 8) (Tq (8 (To (8 (t^ . (4d) 

These vectors characterize the local properties of individual spins, whereas the other tensors describe the spin corre- 
lations. All correlation functions are in the limits 

- 1 < Rc^p^s < 1. (5) 

As idtp"" = [H, p"] (n = 1, 2, 3, . . . ) at unitary evolution there is the numerable number of conservation laws Trp = 
Ci = 1, Tr p'^ = C2, . . . , where C„ are the constants of motion, from which only the first (72,6*3,..., Cie are algebraically 

independent From the conservation of purity, for which (p^)ife = (p)ifc, the polynomial (square-law) invariants 

are obtained. The square polynomials also control the signs RajS-^s- The length of the generalized Bloch vector 
j^epnu conserved under unitary evolution: 

= (6) 

The Liouville-von Neumann equation accepts the real form in terms of the functions Rap-^s a-s closed system of 255 
differential equations for the local Bloch vectors and spin correlation functions 

dtRqOOO = £isqh'^RsOOO + £stq{J'^^RtsOO + J'^"i?tOsO + J'^^Rtoos), (7a) 
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dtRoqOO = Sisqh^RosOO + £stqiJ^^ RstOO + J^^ RotsQ + J^^Rotos), (7b) 

dtRooqo = Sisqh^Rwso + £stq{J™ Rsoto + J^" Rosto + J""Roots), (7c) 
dtRoooq — £isqh^ Rmos + £ stq{J'^"' Rsoot + J^'^Rosot + J"^Roost), (7d) 
dtRqkoo = £isqh^Rskoa + £iskh^ Rqsoa + J'^^£ksq{Rsom ^ Rosoa) 

+£tsqiJ'^"RsktO + J'^^RskOt) + £tsk{J^"'RqstO + RqsOt) , (7e) 

dtRqOkO — £isqhfRsOkO + £iskh" RqOsO + J'^"£ksq{RsOOO ~ RoOso) 

-\-£tsq{J'^^ RstkO + J'^^RsOkt) + £stk{J^^^RqstO + J^^RqOts), (7f) 

dtRqOOk = £isqhfRsOOk + £iskh'^RqOOs + J*^^ £ksq{RsOm ^ RoOOs) 

+£tsq{J'^^RstOk + J'^^RsMk) + £tsk{J^^RqtOs + J^^RqQts), (7g) 

dtRoqkO — £isqh^RoskO + £iskh^i RoqsO + J^"' £ ksqiRosOO ~ RoOso) 

+£tsq{J'^^ RtskO + -P^Roskt) + £stk{J'^^RsqtO + J^^"Roqts), (7h) 

dtRoqOk = £isqh^RosOk + £iskh^RoqOs + J^" £ ksqiRosOO — RoOOs) 

+£ stq{J'^^ RstOk + J^'^Rotsk) + £stk{J'^"'RsqOt + J^^"Roqst), (7i) 

dtRoOqk = £isqh^Rmsk + £iskh^ RoOqs + J""" £ ksq{RoOsO ~ -RoOOis) 

+£stq{J'^"'RsOtk + J^^Rostk) + £stk{J'^"RsOqt + J^^Rosqt), (7j) 

dtRqklO = £isqh1RsklO + Eiskh^RqslO + £islhfRqksO + J*^^ £ ksqiRsOlO — Roslo) 
+ J^^£lsk{RqsOQ ~ RqOso) + J^^£lsq{RskQO ~ ^Ofeso) 

+ J'^"'£stqRtkls + J^"£tskRqslt + £tslRqkst, (7k) 

dtRqkoi — £isqhfRskai + £iskh'-Rqsai + £isih^RqkOs + J'^^£ksqiRsoai — RosOl) 

+ J^^£lsk{RqsOO ~ RqOOs) + J'^^£lsq{RskOO ~ R{)kOs) 

-\-J'^^^estqRskU + J^^£tskRqstl + J^^£ stlRqkst , (71) 

dtRqOkl = £isqhlRsQkl + £iskh^RqOsl + £islhfRqOks + J"^^ £ ksq{RsOOl ~ RoOsl) 
+ J'^"£lsq{RsokO — Rooks) + £lsk{RqOsO — RqOOs) 

+ J'^PetsqRstkl + J^^£stkRqstl + J'^" £ stlRqskt , (7m) 

dtRoqkl = £isqh^Roskl + £iskhfRoqsl + £islhfRoqks + J^^ £ ksq{RosOl " RoOsl) 
''£lsq{RoskO — Rooks) + J^^£lsk{RoqsO — RoqOs) 
+ J'^^£tsqRtskl + J'^^£tskRtqsl + J'^^ £ stlRsqkt, (7n) 

^tRqklm — £isq^i Rsklm ^ £iskhj^Rqslm, ~t~ £islh^ Rqksm £ismhi Rqkls ^ J £ ksq(^RsOlm Roslm) 
^ £lsq{RskOrn ~ Roksm) + J^^''£lsk{RqsOrn ~ RqOsm) + J'^"' £msq{RsklO ~ Rokls) 

-^■P^£msk(,RqslO ~ RqOls) + J^^£rnsl{RqksO — RqkOs)- (7o) 

The derivation algorithm of the system equation has been presented in ■ The set of equations ([7]) with the 
initial conditions has wide applications, since the magnetic field enters in the form of arbitrary functions. It allows to 
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make numerical calculations for: 

1 . continuous (paramagnetic resonance in a continuous mode) , as well as for pulse modes (nuclear magnetic resonance) ; 

2. by means of this system it is possible to investigate the entanglement dynamics of qubits in a magnetic field [l5| 
as the entanglement measures are expressed in terms of the reduced density matrices or of populations; 

3. the important apphcation of the system ([7]) is quantum approach to the Carnot cycle [l[, [H, Q, [1], 
when a working body is a finite spin chain. 

In case of equivalent particles at = hP = h"^ = from the equations (O it follows, that the square length of the 
total magnetization (i?qooo + Roqm + ^oogo + -Rooog)^ is conserved. 
In the external field the energy of system is defined by formula 

= Tr H{t)p = -{h^Riooo + h^Rmm + h^Roam + h^Romi) 

+2{J'^PRiiQ0 + J'^"RiOiO + J^"RoiiO + J'^^RiOOi + J^^RoiOi + J^'^RoOii)- (8) 

The change of the total energy expectation value is equal to dtE — dfTr H{t)p ~ Tt dtH{t)p + Ti H{t)dtp. Due to 
the equation of motion © we have Tt H{t)dtp = —iTrH{t)[H{t),p] = 0. For the external dc field the hamiltonian is 
independent of time, hence dtE = 0, that is the energy of system is the constant in time. 

The change of work W can be associated with the term where only spectrum changes dtW ~ Tr dtH{t) p. The change 

of heat Q is then dtQ = Ti-H{t)dtp. We have, dtdtW = dtTv dtH\t) p^ (integration by parts)= TrHlt)p{t) - 

Tt H{0)p{0) — dtTr [H\t), dtp] = AWsys- The last integral is equal zero due to the equation of motion. Therefore 
the work production system is given by expression of the form 

AWsys = TTH{t)p{t) - TtH{0)po. (9) 

The change of the full system energy in time t which consists of the performed work and heat energy, in our problem 
it is equal to the work production, as in the closed system heat energy is not produced. 

In system assumin g, fo r example, J*^" ^ = J"" = 0, we get the closed system of equations for the description 
of three-qubit dynamics [la|. Assuming J*^" = J^" = J*^" = JP" = 0, we get the closed system of equations for the 
description of two-qubit dynamics 

dtRqo = eiigh^Rio + emiqJ^^Rim, (10a) 

dtRoq — SHqh^Rw + £mlqJ'^^Rmh (10b) 
dtRqk = eaqh^Rlk + SimkhfRqm + J^^ekmqiRniO - Rom) , (lOc) 

where Rqo = Tt pcjq ® (Tq, Roq = Tt pao ® aq, Rkq = Tt pak ® <Jq- 



IV. QUANTUM THERMODYNAMIC VARIABLES 



The reduced density matrices describe dynamics of subsystems (4 matrices of individual particles, 6 matrices of 
two particles, 4 matrices of three particles) and, for example, for a particle e, for the coupled particles ep, epn can be 
written as 

e rp 1 n 1 ( 1 + ^3000 -RlOOO ^ 4-^2000 \ C11„^ 

p _ Trp„„p - -i?„oooa. ~ 2 1, i?iooo + »i?20oo 1 - i?3000 ) ' ^^^^^ 
pP" = ^Rofjjocrp «) CT^, (lib) 

p^P" = ii?a^^oO-Q ®Ufi® G^. (11c) 
o 

The matrices (|lip are determined by the system solution (O, as the equations for the reduced matrices are not closed. 
From the system ([7]) it follows, that probability of spin flip e from initial state is equal to 
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and for the spin p 



It is known Q that in weakly coupled systems in a longitudinal field h*^ = (0, 0, hf) the transverse coherences 
-RiooO:^20oo are equal to zero for the initial state (jA.l[) . that is the reduced matrix p'^ is diagonal in the course of 
unitary evolution. In this case it is possible to define the local dynamic temperature correctly. The local or dynamic 
temperature for two- level system can be defined according to the basic meaning of a thermal state [l^ , : 

T%t) = — (13) 

Po 

where fl^(t) is the transition frequency in the e two- level system, equal to h^. Further we set the Bolzmann constant 
kg equal to 1, therefore all temperatures in energy units. 

The von Neumann entropy after multiplication on Bolzmann constant can then be taken as the thermodynamic 
entropy. The e-spin entropy is equal to 

S%t) = -Trp'lnp' = -pglnpO -Pilnp?, (14) 
where the local populations are equal to 

Po = 2 ' ^1 " 2 ■ ^ ^ 

The time derivative of the entropy for e-spin is defined by the formula 

dtS =f7 = In— ^ . (16) 

The stated concerns each particle. The calculation of the work which are carried out by subsystems, by means of 
ST diagrames is applied in the the papers it is shown numerically, that the entropy S and temperature T are 
dependent thermodynamic variables. For the closed trajectory in ST plane the change of full energy AWsys is equal 
to zero and consequently the area captured by the closed curve in ST plane, determines the work during a reversible 
cycle 

AW = - TdS = - Jj Tit)ait)dt, (17) 

where tc is the duration of a cycle, and the sign is defined according to the rule, if at path tracing clockwise the area 
is situated on the right it obtains a minus sign (heat pump). The spin system is isolated and consequently, the Carnot 
cycles can only refer to sub-systems of the 4-spin system. We use pT|) for e,p,n,u particles. 

The energy of the coupled particle e in a magnetic field /i§ in an environment of three others is equal to ^/i|i?3ooo- 
We use the formula 9((|/i§i?3ooo) = ^(^t/iDiiaooo + ^h^dtRsooo- The work production (the heat production) by spin 
e during cycle is equal to J*^ dt^{dth^) R3000 ( Jq" dt^hldtRsooo)- During cycle the energy change of the particle e 

J^" dtdti'^h^RsQOo) is equal to zero, hence J q" dt^{dth^) R3Q00 =— J^^ dt^h^dtRsooo- Having inserted the equations 
(HH), dini) in dni) we conclude that 



AW = - j^T'dS' ^ J^' dt^idthDR^ooo- 



Thus the definitions of temperature (|13|) , entropy ()14p and work ()17p are coordinated with the work / heat production 
for parameters, for which the ST plots are closed. It concerns also p, n, u particles. 
The Klein- von Neumann inequality looks like 

m 

-Trplnp < -^Pij lnpii < Inm, (18) 

i=l 

where — X^I^i Pa Pa is the diagonal entropy d, m is the number of system states. For the initial diagonal state the 
diagonal entropy possesses property d(0) < d{t) = — Yl^i Pa ^'^Pu uM- 
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The dynamics of a purity measure P = Tt p"^ is connected with the dynamics of entropy S' = — Trplnp as follows. 
If the entropy is equal to zero, the system is in a pure state. At the maximum entropy the system is in the maximally 
mixed state. The purity P has the maximum value 1 for a pure state and the minimum value in the mixed state, 
equal to where m is the number of accessible states. The subsystem purity is expressed in terms of the square 
length of the local or generalized Bloch vector ([6|) : 
the purity for the subsystems pn, epnu is equal to 

pup = 1(1 + 5"P2^^ pepnu ^ J_(^ ^ ^epn«2^^ ^.^g^ 



accordingly, where 6"^ — yY^^j=i ^oijo length of the generalized Bloch vector of the pn system. 

Let's define the entanglement measure p and n spins according to [Toj on system solutions, having entered the 
two-particle entanglement tensor: 

"moijo = Roijo ^ RoiooRoojo- (20) 

The tensor moijo is equal to zero, when the two-particle correlation function i?oijo is factorized in terms of the local 
Bloch vectors (|4b[) . (|4c|) and thus the matrix will be separable, i.e. p^" = ff ® ff . By means of tensor wc shall define 
a measure of the two-particle entanglement in the pn subsystem 



3 E "^o.,o- (21) 

This measure is equal to zero for separable state and it is equal to 1 for the Greenberger-Horne-Zeilinger maximally 
entangled state. This measure is applicable both for the pure and mixed states (in all 6 two-particle measures). 



pn 
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V. NUMERICAL RESULTS 



The quantum thermodynamic devices are subdivided into heat pumps and heat engines depending on functional 
purpose. In our model absorbed or made by system work arises due to displacement power levels by a magnetic field 
po| , [111 , , and it also depends on an initial state of the system. 

We will consider the influence of the variable magnetic field 

h^{hi,h2,h^) (22) 

on dynamics of system, in a case, when the fields operating on spins e, p, n, u arc equal to h'^ = = h"' = = h, 
where hi = 0, h2 ~ 0, ~ 1 + s'mujt, with frequency of an external field lo = 0.04 and the exchange constants are 
equal to 0.01. 

For the parameters a set-valued initial state (|A.ip ( see the Appendix), let's choose the following values for our 
numerical results: /i^(0) = /if(0) = /i^(0) = ^(0) = 1.5, r'=(0) = TP(0) = T"(0) = 0.2, r"(0) = 0.6. In Fig. [1] 
the parametric dependence on entropy and temperature for a cycle are shown, i.e. ST cycles with negative work 
(bold line) which are identical for e,p, u spins. The spin n produces some work (thin line) from {e,p, u)-environment, 
while each of the spins e,p,u (bold line) absorbs work. The work performed by n spin, is equal to the area limited 
by the closed curve according to the formula ([T7)) is equal to 0.1331514. e,p,u spins perform the work equal to 
3x (-0.0443838) = -0.1331514. Thereby the full energy change of 4 spin system AWsys ® is equal to zero in full 
compliance with the general results for the isolated system |23l-|27l|. It confirms the use of temperature and 
entropy (|14p as the effective thermodynamic characteristics. The feature of this initial state and modulation magnetic 
field is that these cycles are repeated without any deformation. (It is known the system equation ([7]) with the periodic 
coefficients according to the Floquet theory has periodic or quasi-periodic solution i.e. the Floquct theory does not 
exclude also the periodic solutions. It depends on a set of the coefficients. We have presented this set.) In other words 
each cycle comes to the end returning to the same initial state. The closure of the ST plots does not depend of the 
amplitudes of the driving field (the form and the area vary only), but it critically depends on frequency to llOl and the 
module of exchange constants, that does not depend of ferromagnetic or antiferromagnetic working gas [2^. After 
replacement of a frequency sign u! the circulation direction becomes opposite for all particles. If the value of parameter 
is T"(0) > 0.2, then the particle n makes work, and the spins e,p,u absorb work. The purity of all system decreases 
with increase of T"(0), but cycle-after-cycle and periodicity remain as it is described, and the areas characterizing 
work, increase with preservation of the algebraic sum which is equal to zero. For r"(0) < 0.2 the spins e,p, u produce 
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FIG. 1: ST diagrams for the separable mixed initial state. All spins are in an identical field h and the exchange constants 
are equal to O.Of. It is seen that the field h — (0,0, f.5(i + sincdt)) causes periodic cycles with parameters of the initial state 
hsiO) = /ig(0) = /i3(0) = h^O) = i.5, r=(0) = TP(0) = T"(0) = 0.2, r"(0) = 0.6, uj = 0.04. The duration of a cycle is equal 
to tc = 2-k/ui = 157.08. The arrows indicate the direction of circulation. In addition to the direction circulation the signs +, — 
specify the negative, the positive work of particles. The vertical arrows in the plots specify the return points. 

the work , and the spin n absorbs the work. If r" (0) approaches to 0.2 the absorbed and produced works approach 
to zero, as at the initial moment there is no temperature gradient in the system (passive state [29j). The numerical 
calculation also shows, that for the parameters of the initial state T^(0) = TP{0) = 0.2; r"(0) = T"(0) = 0.25 (other 
parameters as in Fig. 1) particles e,p absorb work equal to 2(- 0.002265), and particles n,u produce the same work 
2(+0.002265). 

We would Hke to indicate, that in a vicinity t « 117.8 the local temperature of all particles goes to zero as the 
frequencies n^{t) for all particles go to zero. It is necessary to notice, that for the minimum local temperatures 
the eigenvalues of the Hamiltonian ([1]) come nearer to zero, and with the growth of temperatures the eigenvalues 
become bigger. The transition probability of each particle ((T^ is close to 1 and makes one oscillation per cycle. The 
population pieie is approximately equal to 0.925 during a cycle. For the opposite sign /i|(0) = h^{0) = h^{0) — 
h^{0) = — 1.5, pii « 0.925 the thermodynamic characteristics do not change. In Fig. [1] it is clearly visible, that if the 
initial state of a subsystem is more disordered, a subsystem absorbs/produces more (on the module) work [29j . 

The cycle of each particle is determined by the direct coupling of the given particle with the others and does not 
depend of the coupling constants between other particles. The calculations confirm, that the dependence between the 
energy of particles p, n = ^(^f^oioo + ^"-Roojo) + 2JP"i?oiio ^i^d the entropy S'p^ is also cyclic. This property is 
carried out for any particle pair. 

The additional quantum thermodynamic characteristics of eu and pn subsystems are resulted in Fig. [5] and Fig. |31 
Under unitary evolution the global purity P'^v™ does not depend of time. The purity P^" (|19p has maximum at 
an entropy minimum. The work production of eu subsystem is accompanied by the entropy increase and the purity 
reduction and the inverse process occurs at work absorption. It is visible, that the local entanglement mP"^ (PT|) has 
maximum for a minimum entropy 5^" = — Tr p^" In pP". The entanglement (j21[) between the particles is noticeable in 
the middle of cycle and it is bigger between particles with opposite signs in work, and, as the calculations confirm, 
grows with the disorder increase in system. This entropy is always less or equal to the sum entropies of the individual 
spins in pn subsystem. Also the inequality Klein- von Neumann is carried out (jl8p . i. e. the diagonal entropy is bigger 
or equal to pn subsystem entropy. 
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FIG. 2: Time dependence on purity of all system p"^P"" and on characteristics eu and pn subsystems for one cycle with the 
parameters, corresponding to the ST diagrams in Fig. \T\ 



For the cases presented in Fig. [T]the diagonal entropy o?^" (c?"^") coincides with the sum entropies of the individual 
spins SP + S"" {S^ + 5"). 

The control of calculations was carried out with the help invariants of motion, described in section III, and all the 
correlation functions have been in the limits — 1 < Raj3-yS < 1- 

VI. CONCLUSION 

The closed system equations for the local Bloch vectors and spin correlation functions of four two-level systems 
with the exchange interaction, being in the time-dependent external magnetic field is derived. The invariants of mo- 
tion have been found , which are necessary for the control of computing. The numerical analysis of thermodynamic 
behaviour in 4 spin system depending on the parameters characterizing an initial state and modulation of the driving 
field was performed. It was numerically found, that under unitary dynamics the work production of one part of 
system is compensated by absorption of work produced by the other part. The work production by the subsystem is 
accompanied the entropy growth and vice-versa, the entropy becomes less with work absorption by the subsystem. 
It was revealed, that the ST cycle of each particle is determined by the direct coupling of the given particle with the 
others and weakly depends on coupling constants between the other spins. 

It was numerically shown that the formulas for the temperature ([T^ and the entropy are the effective thermo- 
dynamic characteristics for work calculation in 4 spin system in the course of unitary evolution. 
The research of the "spin gas" properties is necessary for the implementation of quantum thermodynamic cycles in 
spin systems [sol [sH . Taking into account an environment zero balance will be broken and the system will become 
heat engine or heat pump depending on temperatures of heat baths. 
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FIG. 3: Time dependence on the local entanglement msM for subsystems pn (thin line) and eu (bold line) for one cycle. The 
parameters correspond to the ST diagrams in Fig. [1] 



Appendix 

The concrete calculations were executed for the separable mixed initial state Sep-Mix: 

PSep^Mix{0)=p''(E)pP(E)p"(E)p", (A.l) 



where 



2cosh/i|/2T»(0) V exp(/iy2T^(0)) 



,i ^ {e,p,n,u), (A.2) 



T'(0) is the initial temperature, h^^{0) is the initial field. 

The initial state psep-Mix(fi) corresponds to a choice of initial values i?oooo = 1, ^ooo3 = 
— tanh(5; i?oo30 = — tanh7, i?oo33 = tanh7tanh(5, i?0300 = — tanh/3; i?0303 = tanh/3tanh(5; i?o33o = 
tanh/3tanh7; i?0333 = tanh/3tanh7tanh(5; i?3ooo = — tanha; i?3oo3 = tanhatanhJ; i?303o = tanhQ;tanh7; i?3033 = 
tanhQ!tanh7tanh(5; i?33oo = tanhatanh/?; i?3303 = tanhQ!tanh/3tanh(5; i?333o = tanhQ;tanh/3tanh7; i?3333 = 
tanhatanh/3tanh7tanh(5, in which a = /i|(0)/2T'=(0), (3 = /if (0)/2rP(0), 7 = /i^(0)/2r"(0), S = mO)/2T''{0), 
and the other initial values are equal to zero. 
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